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ABSTRACT 

The maximum cosmic ray energy achievable by acceleration by a relativis- 
ts blast wave is derived. It is shown that forward shocks from long GRB in 
the interstellar medium are powerful enough to produce the Galactic cosmic-ray 
component just below t he ankle at 4 x 10 18 eV, as per an earlier suggestion 



( ILevinson &: Eichlerl Il993l ). It is further argued that, were extragalactic long 
GRBs responsible for the component above the ankle as well, the occasional 
Galactic GRB within the solar circle would contribute more than the observa- 
tional limits on the outward flux from the solar circle, unless an avoidance sce- 
nario, such as intermittency and/or beaming, allows the present-day, local flux 
to be less than 10 -3 of the average. Difficulties with these avoidance scenarios 
are noted. 

Subject headings: Galaxy, cosmic rays, gamma-ray bursts 



1. Introduction 



The ultrahigh-energy (UHE) range of the cosmic-ray (CR) spectrum is generally broken 
into three parts: 1) the steeper knee-to-ankle segment (~ 10 15 ' 5 to 10 18 6 eV), 2) the flatter 
"trans-ankle" CR below the GZK cutoff (10 19 ' 6 eV), 3) and trans-GZK CR. Trans-ankle CR 
are probably extragalactic, showing little anisotropy at E < 10 19 ' 6 eV, and some anisotropy 
at E > 10 19 6 eV towards the local supercluster. 

In this paper, we consider the hypothesis of sub-ankle UH ECR origin from long GRB 



( ILevinson fc Eichlerl Il993t IWick et al. 



2004 ; ICalvez et all 120101 ). We show that (a) Galactic 
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7-ray bursts (GRB) are sufficiently good accelerators and sufficiently powerful to account 
for sub-ankle UHECR, but that (b) the UHECR near-isotropy limits the current Galactic 
UHECR output per unit star-forming mass to a value far less than what is energetically 
required to account for trans-ankle extragalactic UHECR by extragalactic GRB. Conclusion 
(b) can be generalized to any hypothetical UHECR source whose rate density, like that of 
long GRB, is in proportion to star formation. This challenges any model in which such 
sources account for all UHECR. 

Many past authors have proposed that GRB make the CR above 10 19 eV, but for this 
ene rgy range there remains the alternative hypot hesis that they come from active galax- 



ies ( IThe Pierre AUGER Collaboration et al.l 120101 ) . Doubts remain that GRB could supply 



the h ighest-energy cosmi c rays. The problems include disparity in the total energetics of 



each (lEichler et al.ll2010l . Part I, and references therein), adiabatic losses, which lower the 
maximum energy should the acceleration be in a compact region, and the isotropy problem. 
The isotropy problem, discussed here and in an accompanying paper (Part III), is basically 
that stars in the Milky Way are distributed anisotropically relative to the Earth. If they 
- or sources similarly distributed - were responsible for UHECR, the UHECR should also 
show anisotropy, the Galactic magnetic field notwithstanding. Part III studies in detail the 
propagation of CR from Galactic GRB and specifically compare with data the expected 
anisotropy, composition, and intermittency behavior. 



2. Particle Acceleration in Outflows: A General Discussion 

To be efficiently shock accelerated, a CR that has crossed the shock toward the upstream 
must be overtaken again by the shock of the order of n = (3/(3$ times, where (3, the velocity of 
the particle in units of c, exceeds (3s, the velocity of the shock. If the shock moves at Lorentz 
factor r Si then a CR that is overtaken by a spherical blast wave of age T at radius Rs ~ PscT 
must have been deflected (by gyration or scattering) through an angle AO > l/Tg, while 
residing upstream, within a time At ~ Rs/ Pen = fisRsI '/3 2 cQ In other words, a necessary 
condition for efficient shock acceleration is that AO / At = ZeB/^/mc > (3 2 c/((3sRs^s)- 
Defining the maximum kinetic energy E max for convenience to be ~ (3^ nax ^ max mc 2 , we may 



^^Here it is assumed that, whether the propagation is stochastic or scatter-free, each reorientation or 
reversal of direction must happen within a CR path length of order (3sRs/(3, as the cumulative time a 
CR spends within a gyroradius r g of the shock is only ~ Rs/ (3c. This assumption may be controversial 
for scatter-free propagation in a perpendicular magnetic field, because there is no rigorous proof to our 
knowledge of impossibility of more prolonged trapping, but neither are we aware of any counterexample in 
view of systematic drift (see below) . 
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write 



E max <ZeB(3 s R s r s (1) 



which generalizes p r evious results for diffusive shock acceleration (DSA) Eichlerl ( I198ll ). 



Forman and Druryl (119831 ) and references therein], and for scatter- free shock drift, when 
rs ~ l|§ Note that scatter-free gyration even in perpendicular shocks, though it gives a 
much thinner precursor than stochastic propagation, cannot in general confine a CR particle 
to a subrelativistic blast wave (e.g. a supernova remnant) in all three dimensions, for the 
particle would generally drift off to the side within a time R s /(3c after gaining the potential 
difference ZeB(3$Rs in energy. Also note that we have neglected adiabatic losses. 

Random scattering, for a given magnetic-field amplitude, changes the CR's direction 
more slowly than undisturbed gyration and, for relativistic shocks, usually makes it harder 
for the shock to catch up with a particle. Therefore if turbulent magnetic field amplification 
(MFA) increases the field strength to a value B rms , equation should not be used with 
B = B rms if the coherent scattering angle is less than 1/Tg, i.e. if the coherence length of 
the field I is less than r g /Ts- 

For random small deflections of 58, the mean free path A is about (3c/Dgg, where the 
angular diffusion coefficient, Dgg, is given by Dgg = (59) 2 /5t ~ [ZeB rms / Pjmc 2 } 2 ^ / (1/ (3c) = 
r~ 2 l(3c, where 5t ~ 1/ (3c is the scattering coherence time over which the particle scatters by 
an angle 59, I is the coherence length of the magnetic field, and r g , in a turbulently enhanced 
magnetic field, is defined as r g = Pjmc 2 /ZeB rms . The condition for efficient acceleration is 
now AO 2 = D ge At > 1/T 2 , or r 2 /l ~ A < (3 s R s T 2 s /(3. Finally, we have 



E max < Ze B rms [IRsPsP] 1 T 5 (2) 
which, for a given field strength, is less than the previous expression for E max when / < r g /Ts- 

This expression for E max implies that, if / < r g /Ts, MFA raises E max only if it raises 
the value of B 2 ms l. Simply tangling the field on a small scale so that its strength varies as 
r] < 1/2, does not raise E max . 

For a self-similar, energy-conserving relativistic blast wave in the interstellar medium 



2 Note that the expression E max = ZeB Rs, often taken from figure 1 of Hillas (1984), is consistent with 
equation ((T|) only if f3 s and T s are both of order unity. 
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flBlandford fc McKeeHl976[ ) 



R s « (lT^/ieTrpc 2 ) 1 / 3 ^ 273 « (6 x 10 18 cm) f ^) 1/3 r~ 2/3 . (3) 

where the total energy of the blast is 10 54 -Es4 ergs and the ambient nucleon density is 
nocm -3 . This suggests, again taking the limit for AO/ At as arising from coherent gyration, 
that the gyroradius of a maximally energetic escaped particle is 

(E \ 1/3 

r 9 ,m ax « T S R S ~(6x 10 18 cm) l-^j . (4) 

In the early stages of a powerful GRB blast wave, 0.1 < E 5i < 10, T$ ~ 10 3 , while 
10 -2 < no < 1. Escaping particles, therefore, should obey 10 18 ' 5 < r g < 10 20 ' 8 cm, and 
contribute to the Galactic CR component in the corresponding energy range. They should 
be represented in the flux we observe and in the quantity wq, which is defined below to be 
the CR power per unit baryon mass within the Galactic solar circle. 

The range 10 18 ' 5 cm < r g < 10 20 ' 8 cm in the Galaxy corresponds to the energy range 
10 16 Z[B/10fiG] < E < 10 18 - 3 Z[B/10/iG] eV, precisely the range, to within the uncer- 
tainties, of the "knee-to-ankle" portion of the CR spectrum, which is said to evade the 
capabilities of supernovae. Relativistic blast waves in the Galaxy fill in this range nicely. 

Ultrarelativistic shocks are likely to be quasiperpendicular in the frame of the shock, 
and, on these grounds, their ability to accelerate particles efficiently has been questioned. We 
agree that it is, a priori, a fair concern, but note that nonthermal spectra in GRB afterglows 
seem to indicate that shock acceleration works there just fine. The many ways to evade the 
arguments against shock acceleration in quasiperpendicular geometries are not the subject 
of this paper. 



3. The CR power per unit baryon mass 



If we were to suppose that the mechanism supplying the sub-ankle CR somehow extends 
well beyond the ankle, we would encounter the problem that the Galactic component of these 
CR would be highly anisotropic, assuming their source distribution would be concentrated 
inside the solar circle in the Galaxy, because their transport is no longer fully diffusive and 
includes many Levy flights. This, in addition to the sharp change in spectral index at the 
ankle, is reason to suppose that Galactic GRB limit their output to CR below the ankle. 
Even in the sub-ankle range, the observ ational limits on aniso tropy pose strong constraints 
on the models. Below, and in (Part III: iPohl fc Eichlerll201ll ). this is further quantified. 
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Let f sc ,bM sc = (2 x 10 44 g) f sc ,b be the total baryonic mass within the Solar circle, 
where M sc ~ 2 x 10 44 g. Using the allsky UHECR integral flux F[Ei,E 2 ] in energy interval 
[Ei,E 2 ] (in units of EeV) implied by Auger, and assuming that the fluxes above the ankle 
are extragalactic and uniform in the cosmos, we find that the UHECR source power per unit 
baryon mass in the [4,40] range, u>[4,4o], is 



-^[4,40] - _j _ ls /A[ 440 ]\ 

A[4,4o] MbPc V Gpc / 



where -F[4,4o] = 0.017 erg/cm /yr (lEichler et al.ll2010l ; lAbraham et aU2010al ). and the implied 



luminosity from within the Galaxy's solar sphere, £7[4,4o]> is 

£[ 4 ,40] = ^[4,40],g M sc / sc>6 « (2.5 x 10 38 erg s" 1 ) (^^j /f 1 fsc,b, (6) 

Here A^^y is the "horizon" rangdl of UHECR in the [Ei,E 2 ] range, (A [4,40] ~ 1 Gpc), 
&Bpc ~ 1-4 x 10~ 31 g/cm 3 is the cosmic density in baryons, and f s = w^^/w^q^g] is the 
ratio of the average UHECR source power per unit baryon mass to that in our Galaxy. To 
be precise, the subscript "G' denotes the Galactic value within the solar sphere. Because 
spiral galaxies like our own comprise about half the cosmic mass, with the other half in 
galaxies with less star formation and hence probably lower UHECR source power, one may 
estimate the value of f s to be about 1/2 if UHECR sources are distributed in proportion to 
star formation. 

On the other hand, if the solar system fairly samples the outward flux of cosmic rays 
within the energy range [Ei,E 2 ] through the solar sphere, i.e., if the local flux equals the 
average over the solar sphere, then the inferred power is given by 

E [EuE2] = 4vr F [EuE2] R 2 SC fi [EuEa] = 4tt i& / EP(E)cf(E)dE (7) 

J Ei 

where R sc = 8 kpc is the radius of the solar circle and ^[ El)E2 ]C is the average ratio of enthalpy 
flux [in the anticenter direction, defined to be fi = cos 9 = 1], J'dE J dfi[Efif(E,n)] [where 
the integral runs from E\ to E 2 ], to energy density 4vr f'dE[Ef(E)/c] =2n ft dE J dfi[Ef(E,fi)/c] 
It is measured directly for each energy bin with CR anisotropy measurements. The current 
experimental limits on j3 set by the Auger Observatory are, to 99% confidence, (3 <0.004 in 



3 The horizon range is shorter than the instantaneous range because the expansion of the universe enhances 
the losses both by adiabatic deceleration of the particles and a raising of the background photon energy 
density and the losses it causes in the past relative to the present. 
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the [0.4,4] EeV range and (3 < 0.025 in the [4,40] EeV range (Abreu et al., 2011). Here we 
have used the facts that most of the energy flux is towards the low end of these ranges, where 
the limits on anisotropy are strongest, and that (3 is 1/3 of the first harmonic amplitude given 
by Abreu et al. (2011). Under this assumption, we obtain 



£ [4) 4o] = ^[4,40] 47r J R s 2 c /3 [4i 40] < 1 x 1Q 35 erg s 1 



and correspondingly 



£[0.4,4] = £[o.4,4] 4tt R 2 sc (3 [oaa] < 2 x 10 35 erg s" 1 (9) 
Equations (jSJ) and (jSJ), together with the constraints on (3 imply that 

/^^^^ 2500/^f^iV 1 (10) 

W[4,40],G V G P C / 

Note that all UHECR sources that have a power scaling with star-forming mass, e.g. the 
hypernova scenario for long GRB, should have a high likelihood of being present in the 
Galaxy, i.e. f s < 1. We conclude that if a) the sources of UHECR are fairly represented 
in our own Galaxy, and b) the solar-system location fairly samples these CR at present, 
then the hypothesis that such sources in other galaxies maintain an extragalactic flux at the 
observed level would be inconsistent with the observed CR flux. There would be more CR 
production within the solar circle than allowed by observation. This is a challenge to any 
theory of UHECR origin from long GRB. 



Discussion 



The limit on inferred source power per unit baryon mass required to sustain Galactic 
UHECR in the [4-40] EeV range that is imposed by the observed anisotropy limits is smaller, 
by more tha n 3 orders of magnit ude, than what is required for an extragalactic origin, as 
calculated in lEichler et al.l (120 lOf ). and it corresponds far better to the power per unit mass 
of gamma rays from GRB. This numerical coincidence fits the hypothesis of a GRB origin for 
the Galactic component of UHECR, without invoking a much larger unseen energy reservoir 
for GRB. In fact, it would allow a Galactic origin for UHECR above the ankle were it 
somehow possible to trap these CR within the Galaxy effectively enough to obey the isotropy 
constraint. It remains to be shown that applying the hypothesis of UHECR from Galactic 
GR B to subankle Galacti c CRobeys the isotropy constraint, and this analysis is done in Part 
III JPohl & Eichlerlhoilh . 
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Although the discussion, for historical reasons, has used GRB as a standard for power 
production, it is independent of GRB. The highest-energy CRs, whatever their source, are 
surely extragalactic, and apparently produced with a higher power per unit (star forming) 
mass than that contributed by the matter within the Galactic solar sphere, given the observed 
limits of UHECR outflow from this sphere. This challenges any theory of their origin from 
matter and phenomena of the sort to be found within 10 kpc or so of the Galactic center. 

We have considered several alternative possibilities. AGN are an obvious possibility, as 
they are not represented by our Galaxy, i.e. w/wq ^> 1- 

Another possibility is that the sources are white dwarfs or neutron-star mergers from 
binaries in a very extended halo, and that they spend very little of their time within the solar 
sphere. Conceivably this could include short GRB, although their total energy output in the 
cosmos is probably an order of magnitude less than even that of long GRB, so the question 
of total energetics would still loom large. On the other hand, short GRB, not being tied to 
the SFR, need not suffer the recent decline in rate relative to earlier epochs, and so could 
be an order of magnitude more common, relative to long bursts, at present than in earlier 
epochs. In any case, one would still have to check that the implied flux is below observed 
levels and of a suitable angular distribution. Why, for example, would there be correlation 
above the GZK cutoff with the local supercluster? If one is willing to attribute sub-GZK CR 
above the ankle to a different class of sources from those above the GZK cutoff, then short 
GRB in the Galactic halo may account for the former, provided they are distant enough to 
respect the strong limits on anisotropy. 

In an effort to accommodate the hypothesis of a GRB origin for all UHECR, we have 
also considered the possibility that our present location does not fairly sample the UHECR 
exiting the solar sphere, and that the large UHECR output that would be necessary to 
supply all of the UHECR at energies where their flux would be extragalactic could then 
mostly evade the solar system. They conceivably could, for example, be blown out in jets 
that have avoided our location and/or with an intermittency that excluded the present epoch 
receiving a fair representation of the time average. If GRB blasts were to escort all their 
CR safely out of the Galaxy in narrow jets that avoid our location, there would be less of 
an anisotropy problem associated with a GRB origin for extragalactic UHECR. But this 
scenario would differ from the common view that GRB blasts slow down to subrelativistic 
Lorentz factors, spreading in angle, within the Galaxy. If the UHECR escape the jet, then 
according to equation 4, they probably get significantly deflected before escaping the Galaxy 
at large, and it is not obvious that they could remain sufficiently collimated to conform 
to an avoidance scenario. A single jet, if it leaks UHECR, contaminates the sky with a 
strongly anisotropic component at the energy range in which CR are strongly scattered by 
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the Galactic magnetic field, but not contained by the jet. If, on the other hand, the jet 
contains all the UHECR, then the latter suffer enormous adiabatic losses. The question 
is whether all but ~ of the UHECR can avoid leaking or escaping into the Galaxy 
at large and mixing with its CR population. This would appear to require a scenario in 
which CR at ~ 10 19 eV would be extremely well confined to the shock (strong scattering) 
without suffering adiabatic losses and without being scattered out of the shock's path into 
the interstellar medium. 

Intermittency may explain a low flux from Galactic long GRB, if the time between GRB 
per Milky- Way-type galaxy were more than R sc /(3cj3) ~ 10 4 //3 years, the escape time of CR 
from the Galaxy. The collimation of GRB jets to within several degrees however, which is 
now believed to be the case, suggests that GRB are as frequent as 10~ 7 f^ 1 per year per M sc , 
where /& is the beaming factor, believed to be of order 10 -L5 to 10 -3 . Specifically, the local 
rate of GRB that we detect is R = R\ Gpc _3 yr _1 , with R\ ~ 1. The expected rate Rq 
within the Galaxy should then be 

-9 rrW-l'N fsc,bRl (n - 1n _5 -l\ D f ( fb \ ( fs^B 



R G ~ (1.6 ■ 10- 9 yr- 1 ) J f c > ^ ~ (2.5 • 1(T 5 yr- 1 ) R 1 f sc b I il^L ) . (11) 

G v y hfs^B V 10 v V 0.02 y v ; 

Given the near isotropy of UHECR, if a good fraction of them were produced within the 
solar circle, their escape time from the Galaxy, R sc //3c would be larger than 1/Rg, and 
there would necessarily be many Galactic GRB per escape time. Moreover, if the escape is 
exponential, more that 10 escape times would be necessary to clear all but 10~ 3 of the CR 
released by the GRB. These considerations suggest that supply intermittency from GRB is 
most plausible when j3 ~ 1, but it is doubtful that the assumption of j3 ~ 1 is consistent 
with the strong limits on anisotropy. 



On the other hand, the original scenario of iLevinson fc Eichlerl (119931 ). in which Galactic 



GRB supply the Galactic CR component below the ankle, need not make significant energy 
demands on the GRB, because there is no constraint imposed on the required output per 
unit mass other than what the local Galactic flux dictates. Given the low anisotropy of 
UHECR, it is likely that those below the ankle are confined to the Galaxy for nearly 10 6 
years, in which case the rate of GRB and their energy output are consistent with a CR 
production per GRB that is less than that of gamma radiation per GRB. This however, 
woul d be incompatible w ith the hypothesis that GRB produce UHECR above the ankle as 



well f lEichler et al.ll2010h 



While the low observed anisotropy eases the energy demands on sources of Galactic 
UHECR, it imposes a strong constraint of its own. It remains to be shown that sources of 
UHECR, if distributed as luminous stars in our Galaxy, indeed satisfy the constraint of low 
anisotropy. This question depends on the question of the mean free paths of UHECR, hence 
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on their composition, and is taken up in a companion paper (IPohl fc Eichlerl 1201 ll ) . 

Here we have concluded that GRB are energetically sufficient to provide the sub-ankle 
Galactic CRs, given what is known about shock acceleration, relativistic blast waves, and 
GRB parameters. If the hypothesis of GRB origin for sub-ankle UHECR is true, it may have 
implications for the Galactic magnetic field and/or the distribution of those GRB. 

This research was supported by the Israel-US Binational Science Foundation, the Israeli 
Academy of Science, and The Joan and Robert Arnow Chair of Theoretical Astrophysics. 
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